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Abstract. - The quasienergy spectrum of a Bloch electron affected by dc-ac fields is known 
to have a fractal structure as function of the so-called electric matching ratio, which is the ratio 
of the ac field frequency and the Bloch frequency. This paper studies a manifestation of the 
fractal nature of the spectrum in the system "atom in a standing laser wave" , which is a quantum 
optical realization of a Bloch electron. It is shown that for an appropriate choice of the system 
parameters the atomic survival probability (a quantity measured in laboratory experiments) also 
develops a fractal structure as a function of the electric matching ratio. Numerical simulations 
under classically chaotic scattering conditions show good agreement with theoretical predictions 
based on random matrix theory. 



1. - In this letter we study the spectral and dynamical properties of a Bloch particle 
affected by static and time-periodic forces: 

H = p 2 /2 + cos a; + Fx + F^ cos(wi) x , p — ~h'd/dx , (1) 

where Ti is the scaled Planck constant (see below) . Originally this problem was formulated for a 
Bloch electron in dc-ac electric fields H — p 2 /2m + V{x) + e[E+E u cos(wi)] x, V(x+a) = V(x) 

and attracted much attention because of the similarity with the Hofstadter problem 
Indeed, the energy spectrum of a Bloch electron in a 2D lattice under the action of a constant 
magnetic field B depends on the magnetic matching ratio f3 = h/eBa 2 (a is the lattice constant) 
and has a fractal structure as function of this parameter. Analogously, the quasienergy 
spectrum of system ([[]) depends on the electric matching ratio 

u 2irF 
a = — > w B = —j- , (2) 
n 



where ujg is the Bloch frequency, which is u>b = eEa/h in the case of a crystal electron. For 
rational ratios a = r/q the quasienergy spectrum has a band structure; it is discrete, however, 
for irrational values of a Q. Numerically the quasienergy spectrum of the system (Q) was 
studied in ref. Q by using the tight-binding approximation. It was show that for a — r/q 
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the quasienergy bands are arranged in a structure resembling the famous Hofstadter butterfly 
(see Fig. 1 in Ref. ||). It should be pointed out, however, that the results of paper Q only 
partially describe the spectrum of the system (Q) because the tight-binding approximation 
neglects the decay of the quasienergy states. The actual quasienergy spectrum is complex, 
where the imaginary part of the spectrum defines the lifetimes of the metastable quasienergy 
states. 

Because of the extremely small lattice period in crystals, the fractal structure of the 
quasienergy spectrum has never been observed in solid state systems. However, a signature 
of it was recently found in an experiment with cold atoms in an optical lattice Q. The latter 
system models the solid state Hamiltonian (Q), where the neutral atoms moving in the optical 
potential V(x) (Ul is the laser wave vector) take over the role of the crystal electrons. The effect 
of the electric fields can be mimicked, for example, by the inertial force induced by accelerating 
the experimental setup as a whole. (In practice, however, the acceleration was obtained by an 
appropriate chirping of the laser frequency.) The system "atom in a standing wave" has an 
essentially larger lattice period than the solid state system and is, in addition, free of relaxation 
processes due to scattering by impurities and the Coulomb interaction. These features of the 
system were utilized earlier in ref. j^] to observe experimentally the Wannier-Stark ladder of 
resonances. The main modification of the experiment H in comparison with the experiment 
[[| is that a strong periodic driving with frequency to = (r/q)ujB was imposed. (Only the cases 
r/q = 1/2 and r/q = 1/3 were reported.) Then the atomic survival probability as a function 
of the frequency of the probe signal shows additional anti-peaks (see fig. 3 in ref. Q ) , which 
were interpreted as an indication of the fractal structure of the spectrum. We note that the 
width of these anti-peaks is given by the width (i.e. the inverse lifetime) of the first excited 
Wannier-Stark resonances. This imposes a fundamental restriction on the resolution and it 
seems impossible to resolve matching ratios a = lo/lob for a different from lowest rational 
numbers. 

2. - The discussed papers H |J study the system in the deep quantum region. In the present 
paper we discuss the manifestation of the fractal structure of the quasienergy spectrum in the 
semiclassical region of the system parameters. We shall show that in this region the fractal 
nature of the spectrum can be observed without using a probe signal and the electric matching 
ratio can be resolved with any desired accuracy. 

The characteristic measure of the systems "classicality" is the scaled Planck constant ft' 
entering the momentum operator in the Hamiltonian ([!]) . Referring to the system "atom in a 
standing wave" the scaled Planck constant is given by 




where uj rec — ftk 2 /2m is the atomic recoil frequency and Vq is the depth of the optical potential 
V(x) — Vq cos 2 (fc£ir). In the experiments Q and || the value of the scaled Planck constant 
was ft! ps 1.5 and 1.6, respectively. In our numerical studies we use ft! = 0.25. Since the 
amplitude Vo of the optical potential is proportional to the square of the laser field amplitude, 
this implies a larger intensity of the laser. 

We simulate the wave packet dynamics of the system (Q) by numerical solution of the time- 
dependent Schrodinger equation in the momentum representation using parameters ft' = 0.25, 
uj = 10/6, F u = 4.16 with values of the static field F in the interval 0.030 < F < 0.083. Then 
the survival probability 

P(t)= [ |Vb,i)| 2 dp, (4) 

J\p\<Po 
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Fig. 1. - Time evolution of the atomic survival probability P(t) calculated for rational values of the 
electric matching ratio uj/uib = J"/g, ? < 420 (see text). With increasing time t = 40 Tb,. ■ -t — 160 Tb 
the survival probability develops pronounced peaks at rationals with small denominators. The system 
parameters are ti = 0.25, u) = 10/6, F^/uo 2 = 3/2, and 0.031 < F < 0.077. 



is calculated. The initial state is the localized Wannier state, which for the cosine potential 
practically coincides with a minimal uncertainty wave packet centered at x = 7r. The moti- 
vation for the chosen value of po — 6, which is much larger than the region of support of the 
initial wave packet, is given in sec. 3 below. This numerical simulation models the experimental 
situation of ref. with the modification that there is no probe signal but the amplitude F of 
the static force is varied. 

The result for the survival probability as a function of F — or, equivalently, as a function of 
the matching ratio a — uj/lob — is shown in fig. ^ and ^(a), which is the central result of the 
paper. The dots in fig. ^ mark the survival probability P{t) calculated for increasing values 
of the observation time t = 40 Tb,. . .t — 160 Tb- The values of static force F are chosen 
for rational values of the electric matching ratios lo/ujb = r/q with q < 420. (Explicitly, 
we considered for q all divisors of 2 2 • 3 • 5 • 7 = 420 and all coprime r- values in the interval 
6/7 < r/q < 15/7.) To guide the eyes, the dots arc connected by a solid line. The actual width 
of the peaks, which is inversely proportional to the observation time, is smaller in fig. ^|(a), 
which shows a more fully developed distribution for a longer time t — 200 Tb- It is seen in 
fig. |^(a) that there are pronounced peaks above an almost constant background which appear 
at low order resonances between the driving frequency and the Bloch frequency. The seven 
largest peaks in the figure correspond (from left to right) to r/q = 1, 5/4, 4/3, 3/2, 5/3, 7/4, 
and r/q — 2, and the observed peak-heights fall off rapidly with increasing denominator q. As 
illustrated by fig. [I], this peak structure develops gradually in time, starting from P(0) = 1 at 
time zero and originates from the different long-time behavior of the survival probabilities. We 



4 



EUROPHYSICS LETTERS 



find exponential decay in time for irrational values of uj/ujb and algebraic decay for rational 
ones [||. Thus the survival probability as a function of F reflects the fractal nature of the 
spectrum in the long-time regime. Moreover, there is no fundamental resolution constraint 
and an arbitrary number of peaks can be resolved by increasing the observation time. We 
also note that instead of varying the static force one can vary the frequency w of the driving 
force. In this case, having in mind a laboratory experiment, it looks reasonable to use the 
gravitational force as the static force JtJ . 

3. - The key point of the conducted numerical experiment is that the amplitude F u and the 
frequency lo of the driving force are chosen such as to insure chaotic dynamics of the system in 
the classical limit. We furthermore choose po in (0) to be larger than the boundary between 
the regular and the chaotic component of the classical phase space (see fig. 1(b) in ref. ||). 
Then the classical survival probability decays exponentially H 

P cl (t) = exp(-ut) , (5) 

where the decay coefficient is, the inverse classical lifetime, is determined by the classical 
Lyapunov exponent and the fractal dimension of the chaotic repellor. We computed the 
classical decay rate numerically and determined the F-dependence of the classical decay rate 
as v sa 0.15 F in the parameter region considered here. It should be noted that the exponential 
decay of the classical probability cannot be considered as an universal law. In some systems it 
is a transient phenomenon and changes to an algebraic decay caused by long-lived trajectories 
sticking near stability islands ^ p0[ . However, this is not the case for system ([!]) and no sign 
of an algebraic decay of the classical probability was detected (at least until time t — 200Tb, 
which was the maximal time in our numerical simulation). 

Because we measure the observation time t in units of the Bloch period Tb = U /F, the F 
dependence in Eq. (||) cancels and the classical survival probability is practically constant. The 
quantum results follow closely the classical exponential decay (||) for irrational values of the 
matching ration lo/lob providing the flat background of the quantum results shown in fig. |[ 
The peaks above the classical plateau for resonant driving, i. e. rational values of the matching 
ratio a = lo/lob are a quantum "stabilization" -phenomenon, which can be understood as 
follows. 

It was shown in ref. JllJ that the eigenvalue problem for the quasienergies of the system (Q) 
for a = r/q can be mapped onto an effective scattering problem with q open channels. When 
the matching ratio a is an irrational number, the number of channel is infinite and the system 
follows the classical dynamics provided the condition %' -C 1 is satisfied. When the matching 
ratio is a rational number, however, the number of decay channels is finite and the quantum 
system (independent of the value of ft!) is essentially more stable than the classical one 
In this case the behavior of the survival probability P(t) differs from the exponential decay 
(^|) and is determined by the distribution of the imaginary parts of the quasienergies, i.e. the 
distribution of the resonance widths. 

In the case of chaotic classical dynamics an analytic expression for the resonance statistics 
is supplied by (non-hermitian) random matrix theory (RMT) |l^, fl3| |. The validity of RMT 
for system ([!]) was checked numerically in ref. JfiJ |(| and a satisfactory correspondence was 
noticed. Converting the result of RMT from energy to time domain shows that the decay of 
the probability follows asymptotically an inverse power law ]l4| ] 

P(t) ra {Twt/qtiy , t > t* ra qti/Tw , (6) 

where T\y is the Weisskopf width, which is a free parameter in the abstract RMT. Identifying 
the parameter ITV/Ti' with the classical decay coefficient v, eq. (||) and eq. (||) can be combined 
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Fig. 2. - (a) Atomic survival probability P(t) as shown in fig. gj, however for a longer time t = 200 Tb, 
where the fractional structure is much more clearly developed. To guide the eyes the dots are connected 
by a solid line but the actual width of the peaks, which is inversely proportional to the observation 
time, is smaller, (b) Survival probability P(i) = (1 + vt/q)~ q (eq. (mj) derived from random matrix 
theory for t = 200 Tb- The value of the classical decay rate is v = 0.15 F. 



in the single equation 

P(t)=(l + |) " , (7) 

which has the correct short- and long-time asymptotic and provides a first crude approximation 
to the more elaborate RMT result |1J]. 

Figure §(b) shows the values of the function (0) for t = 200 Tb and the same values of 
the matching ratio a as in fig. ^](a) where we use a slightly different graphic presentation of 
P(t) to stress that the function (m) is a discontinuous function of a for any t. In contrast, 
the atomic survival probability shown in fig. |^(a) is a continuous function of F where its 
fractal structure develops gradually as t — * oo. In fact, the probabilities (^) calculated for 
two close rational numbers u\ and c*2 follow each other during a finite "correspondence" 
time. (For instance, for a.\ = 1 and a 2 — 999/1000 the correspondence time is found to be 
about 50 Tb.) Thus it takes some time to distinguish two close rationals, although they may 
have very different denominators and, therefore, very different asymptotics (^|). With this 
remark reserved, a nice structural (and even semiquantitative) correspondence is noticed. In 
addition, it seems worthwhile to note that also the pronounced quantum resonance peaks, 
i.e. the quantum algebraic decay P(t) « (vt/q)~ q predicted by RMT, is mainly determined 
by the purely classical ^-coefficient due to classically chaotic scattering dynamics. 

4- - We have analyzed the system ([!]) in context with recent experiments studying the 
dynamics of cold atoms in a standing laser wave W. It is shown that in the semiclassical 
region of the system parameters the atomic survival probability as a function of the static 
force (or, alternatively, of the driving frequency) shows a fractal structure. This fractal 
structure is actually related to the fractal nature of the quasienergy spectrum determined by 
the degree of rationality of the electric matching ratio (0). In fact, when a rational sequence 
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of a — r/q converges to some irrational value, the quasienergy bands progressively split into 
sub-bands. This process is accompanied by loss of stability of the quasienergy states and shows 
up, finally, in the complicated (fractal) structure of the atomic survival probabilities which can 
be measured in laboratory experiments. 

Finally, we would like to distinguish the fractal structure of the survival probability discussed 
above from the fractal structure of the survival probability studied in papers [jl5| [if], [jjj . The 
latter appears as a quantum manifestation of the hierarchical island structure of classical phase 
space in a system with an algebraic decay of the classical probability. The origin of the former 
phenomenon, however, is the fluctuating number of the decay channels depending on the value 
of the electric matching ratio (||), which is the control parameter of the system ([!]). 

*** 

This work has been supported by the Deutsche Forschungsgemeinschaft (SPP 470 'Zeitab- 
hdngige Phanomene und Methoden in Quantensystemen der Physik und Chemie'). 
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